Modulational instability in asymmetric coupled wave functions 
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Abstract 

The evolution of the amplitude of two nonlinear ly interacting waves is considered, via a set 
of coupled nonlinear Schrodinger-type equations. The dynamical profile is determined by the 
wave dispersion laws (i.e. the group velocities and the GVD terms) and the nonlinearity and 
coupling coefficients, on which no assumption is made. A generalized dispersion relation is obtained, 
relating the frequency and wave-number of a small perturbation around a coupled monochromatic 
(Stokes') wave solution. Explicitly stability criteria are obtained. The analysis reveals a number 
of possibilities. Two (individually) stable systems may be destabilized due to coupling. Unstable 
systems may, when coupled, present an enhanced instability growth rate, for an extended wave 
number range of values. Distinct unstable wavenumber windows may arise simultaneously. 
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Amplitude modulation (AM) is a widely known nonlinear mechanism dominating wave 
propagation in dispersive media it is related to mechanisms such as modulational insta- 
bility (MI), harmonic generation and energy localization, and possibly leads to soliton for- 
mation. The study of AM generically relies on nonlinear Schrodinger (NLS) type equations 
Q|; a set of coupled NLS (CNLS) equations naturally occurs when interacting modulated 
waves are considered. CNLS equations are encountered in physical contexts as diverse as 
el—gnetic wave propagation QQ, optical ^ DD P la sm a waves BflQ, teaae- 
mission lines ^(J, and left-handed (negative refraction index) metamaterials (LHM) [llj]. A 
similar mathematical model is employed in the mean-field statistical mechanical description 
of boson gases, to model the dynamics of Bose-Einstein condensates 3, Q • In this pa- 
per, we shall investigate the (conditions for the) occurrence of MI in a pair of (asymmetric) 
CNLS equations, from first principles. A set of stability criteria are derived, to be tailor cut 
to a (any) particular problem of coupled wave propagation. 

a. The model. Let us consider two coupled waves propagating in a dispersive nonlinear 
medium. The wave functions (j = 1,2) are modelled by ipj expi(kjr — ujjt) + c.c. (complex 
conjugate), where the carrier wave number kj and frequency Uj of each wave are related by 
a dispersion relation function ujj = ojj(kj). Nonlinearity is manifested via a slow modulation 
of the wave amplitudes, in time and space, say along the The amplitude evolution 

is described by a pair of CNLS Eqs. 

The group velocity v g j and the group-velocity-dispersion (GVD) term Pj corresponding to 
the j— th wave is related to (the slope and the curvature, respectively, of) the dispersion curve 
via v g j = dujj/dk x and Pj = d 2 ujj/2dk^ (differentiation in the direction of modulation) . Qjj 
and Qjj' model carrier self-modulation and wave coupling, respectively. No hypothesis holds, 
a priori, on the sign and/or the magnitude of either of these coefficients. The group velocities 
are often assumed equal, in which case (and only) the corresponding terms are readily 
eliminated via a Galilean transformation. The combined assumption Pi = P 2 , Qn = Q22 
and Q12 = Q21 is often made in nonlinear optics 0, 3|. The case Pi = P 2 , Q11 = Q21 and 
Qio = Q22 was also recently reported, in negative refarction index composite metamaterials 
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b. Modulational (in)stability of single waves. Let us first briefly outline, for later refer- 
ence, the analysis in the case of a single modulated wave, say here recovered by setting ip2 — 



in Eqs. 
ism 



qs ©; t 



0); the (single) NLS equation is thus obtained. According to the standard formal- 

t/>, dropping the index in this paragraph) is modulationally unstable (stable) 

if PQ > (PQ < 0). To see this, one may first check that the NLSE is satisfied by the 

plane wave solution i/j(x, t) = ipo e l Wol 2 *_ The standard (linear) stability analysis then shows 

that a linear perturbation (say ?/>o — > ipo + €o~ipo) with frequency O and wavenumber K [i.e. 

Sip ~ expi(Kx — Qt)} obeys the dispersion relation: (Q — v g K) 2 = P K 2 (P K 2 —2Q \ipo\ 2 ), 

which exhibits a purely growing unstable mode if K < K cr $ = (2Q/P) 1 / 2 \ipo\ (hence only if 

PQ > 0). The growth rate a = ImQ attains a maximum value a max = Q \ipo\ 2 a ^ K cr>0 /y/2. 

For PQ < 0, on the other hand, the wave is stable to external perturbations. 

c. Coupled wave stability analysis. In order to investigate the modulational stability 

profile of a pair of coupled waves, we shall first seek an equilibrium state in the form ipj = 

if)joexp[iipj(t)] (for j = 1,2), where ipjo is a (constant real) amplitude and <Pj(i) is a (real) 

phase, into Eqs. We thus find a monochromatic (fixed-frequency) solution of the 

form fj(t) = Qjot, where Qjo = Qjjipj Q + Qjitpfo (for j ^ I = 1,2, henceforth understood 

everywhere). Considering a small perturbation around equilibrium, we take ipj = (ipjo + 

eif)ji) exp[i<fj(t)], where i[)ji(r,t) is a small (e <C 1) complex amplitude perturbation of the 

wave amplitudes. Substituting into Eqs. (0) and separating real and imaginary parts by 

writing ipj% = a,j + ibj (where dj, bj G 9ft), the first order terms (in e) yield 

db ■ db ■ d 2 a ■ 

- -jfi ~ v 9,3-q^ + P j-q^J + -Qji l -%":, + 2Q i ^ ^oaz = , 

daj daj „ d 2 bj 

Eliminating bj, these equations yield 

kE + V "ai) +^(^^ + 2Q 11 A)^j Ol + P 1 Q 12 ^ 20 ^a 2 = 0, (3) 

(along with a symmetric equation, upon 1 <-> 2). Let Oj = a j0 exp[i(i^x — Ot)] + c.c, 
where K and Q are the wavevector and the frequency of the perturbation, respectively, viz. 
d/dt — > — iQ and — > iK. We thus obtain an eigenvalue problem in the form Ma = Q 2 sl, 
where a = (aio, ct2o) T and the matrix elements are Mjj = PjK 2 (PjK 2 — 2Qjjtp 2 ) = Q 2 and 
M jt = -2P 3 Q^ 3 ^K 2 (for l± j = 1 or 2) 

[(n - <v#) 2 - ng [(n - ^ 2 ^) 2 - nl] = n« (4) 



where = M12M21. This dispersion relation is a 4th order polynomial equation in Q. 
d. Equal group velocities. For v 9j x = v 9t2 , setting fl — v g> i/ 2 K — > Q reduces flU) to 

n 4 - m 2 + D = , (5) 

where T = TrM = Q 2 + Q 2 and -D = DetM = f2 2 f2 2 — Q 4 are the trace and the determinant, 
respectively, of the matrix M. Eq. (J5J) admits the solution Q 2 = |[T ± (T 2 — 4P) 1 / 2 ], or 



^ = + ^2) ± 5 [(«? " ^) 2 + ^ • (6) 

Stability is ensured (for any wavenumber K) if (and only if) both of the (two) solutions 
of (|5|). say are positive (real) numbers. In order for the right-hand side to be real, the 
discriminant quantity A = T 2 — AD = (f2 2 — fi 2 ) 2 + AQ 4 has to be positive. Furthermore, 
recalling that the roots of the polynomial p(x) = x 2 — Tx+D, say r = r 12 , satisfy T = ri + r 2 
and .D = r\r 2 , the stability requirement is tantamount to the following three conditions being 
satisfied simultaneously: T > 0, D > and A = T 2 - AD > 0. 

The first stability condition, namely the positivity of the trace T: T = K 2 [K 2 J2j Pf — 
PjQjjip"jo\ > 0, depends on (the sign of) the quantity q\ = J2j PjQjjl^jol 2 which has 
to be negative for stability. The only case ensuring absolute stability (for any ipj and k) is 

PtQn < and P 2 Q 22 < . (7) 

Otherwise, T becomes negative (and thus either VL 2 _ < < Vt 2 + or Q 2 _ < Q\ < 0) for K 
below a critical value K„± = (2J2j PjQjjtfo/ J2j Pf) 1/2 > (cf. the single wave criterion 
above); this is always possible for a sufficiently large perturbation amplitude \ip2o\ if, say, 
P2Q22 > (even if P\Qn < 0). Therefore, only a pair of two individually stable waves can 
be stable, or the presence of a single unstable wave may de-stabilize its counterpart. 
The second stability condition, namely positivity of the determinant D, amounts to 

D(K 2 ) = P^K^P.K 2 - 2QmP 2 l0 ) (P 2 K 2 - 2Q 22 ^ 2 ) - AQ 12 Q 21 ^ 2 20 ] > . 

We see that D(K 2 ) bears two non-zero roots for K 2 , namely K 2 :)1 i D2 = 2P lp 2 [<?2 T (<?| — 
4P 1 2 P 2 2 q , 3) 1 / 2 ], as obvious from the expanded form D = k 4 (PiP 2 K 4 — q 2 K 2 + q%), where q 2 = 
EKhj = 2PiP 2 (P2Qii^o 2 i + ^iQ22^o 2 2 ) andg 3 = T\K 2 D ^ = 4P 1 P 2 (Q 11 Q 22 -Q 12 Q 21 )V> 2 x ^ 2 . 
The condition D > (V-K" G 3?) requires either: 

- that the discriminant quantity A' = q 2 - AP 2 P 2 q 3 = AP 2 P 2 [(P1Q22VI0 ~ Wii^o) 2 + 
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4P\P 2 Q\ 2 Q ^lV'ioV'io] be non-positive, i.e. A' < (this is only possible if P1P2Q12Q21 < 
and for a specific relation to be satisfied by the perturbation amplitudes ipjo; that is, it 
cannot be generally satisfied, V^o), or 

- that A' > and both of the (real) non-zero roots Kj )1 , D2 of D(K 2 ) be negative; this is 
ensured if q 2 < and g 3 > 0. If q% > and 02 > 0, then the two roots K 2 D l , 2 will be positive 
(0 < if f) ! < if|) 2 ) an d the wave pair will be unstable to a perturbation with intermediate 
if, i.e. iff, ! < if 2 < K 2 D2 . If q 3 < (regardless of g 2 ), then If |, x < < if|) 2 , and the wave 
pair is unstable to a perturbation with if 2 < if|, 2 . 

These instability scenaria and wavenumber thresholds are sufficient for symmetric wave 
systems (i.e. upon 1 <-> 2), as we shall see below. 

The last stability condition regards the positivity of the discriminant quantity A = T 2 — 
4D (irrelevant if D < 0). We consider the inequality 

A(k 2 ) = if 4 (d 4 if 4 - d 2 K 2 + do) > 

where d 4 = (P 2 - P 2 ) 2 , d 2 = 4(P 2 - P^Qn^ - P 2 Q 22 ^ 2 ) and d = A[(P 1 Q 11 ^ 2 l0 - 
P 2 Q 22 ip 20 ) 2 + 4PiP 2 Qi 2 Q 2 i'ij) 2 'ij)2 ]. We should distinguish two cases here. 

If Pi = P 2 = P, this condition reduces to d > 0, i.e. here d' = 4P 2 [(Qnip 2 — Q 22 ifj% ) 2 + 
^QnQ 21^10^0} = 4P 2 g 4 . Stability (for all ipj ) is thus only ensured if Q\ 2 Q 2 \ > 0. For 
symmetric wave pairs, i.e. for P\ = P 2 = P and Q\ 2 = Q 2i , this last necessary condition for 
stability is always fulfilled. If, on the other hand, < 0, the wave pair will be unstable 

in a range of values (e.g. of the ratio ^10/^20)? to be determined by solving d' < 0). 

Let us now assume (with no loss of generality) that P± > P 2 . Since A" = d\ — 4d^do = 
-64PiP 2 Qi2<22i(i : i 2 - Pffii'io'tp'ioi the stability condition A > is satisfied for all if and 
ipjo only if the quantity g 5 = P1P2Q12Q21 is positive, hence A" < 0; again, this is always 
true for symmetric waves. Now if, on the other hand, g 5 < (i.e. A" > 0), then one needs 
to investigate the signs of d 2 = if \ : + if^ 2 = q 6 and d = K\ X K\ 2 = q 7 , in terms of 
the amplitudes ipj . Here, ifi,i/2 = lil^ 2 T \J^2~ 4d 4 d ] . Similar to the analysis of the 
previous condition (see above), one may easily see that both signs are possible for both 
quantities d 2 and do. The only possibility for stability (Vif ) is provided by the combination 
d 2 < and d > (hence if^ ;1 < if^ 2 < 0). The possibility for instability arises either for 
ifi ;1 < < if 2 < ifi 2 (if d < 0), or for < ifj^ < if 2 < K\ 2 (if d > and d 2 > 0). As 
above, see that we obtain the possibility for a window of instability far from if = 0. 
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Instability is manifested as a purely growing mode, when one or more of the above 
conditions are violated. In specific, if T < and/or D < 0, then one (or both) of the 
solutions of the dispersion relation (|1J) (for Q 2 ) becomes negative, say Q 2 < [given by JBJ)]; 



the instability growth rate in this case is given by a = J — fi 2 ., and is manifested in the 
wavenumber ranges [0, -fC cr ,i] and either [0,^^,2] or [Kd,i, Kd,2] (depending on parameter 
values; see the definitions above). 

If A = T 2 — AD < 0, on the other hand (hence D > 0), then all solutions of (@J) are 
complex, thus developing an imaginary part In^f^) = ±v/jA|/2, so (the maximum value 
of) Im(f2±) = Im^^) 1 / 2 gives the instability growth rate a. As found above, this will be 
possible for wave numbers either in [0, K/s.,2] or [i^A.i, -^,2] (see the definitions above). 

The analysis indicates that up to three different unstable wavenumber "windows" may 
appear; these windows may either be partially superposed, or distinct from each other. One 
may therefore qualitatively anticipate MI occurring for K G [0, K cr ] (some threshold) and, 
also, for K e [K' cr , K" r ] [K' cr may be higher or lower than K cr , depending on the problem's 
parameters). Furthermore, the instability growth rate witnessed may be dramatically mod- 
ified by the coupling, both quantitatively (higher rate) and qualitatively (enlarged unstable 
wavenumber region); see in Fig. ^ 

Summarizing the above results, should one wish to investigate the occurrence of modu- 
lational instability in a given physical problem, one has to verify condition ((7j), and then 
consider the (sign of the) quantities qi, ...q 7 (defined above). 

e. The role of the group velocity misfit. It may be interesting to discuss the role of the 
group velocity difference, in a coupled wave system. Keeping the discussion qualitative, we 
shall avoid to burden the presentation with tedious numerical calculations. One may rather 
point out the role of the group velocity misfit via simple geometric arguments. Inspired by 
an idea proposed in Ref. |7|, we may express the general dispersion relation Q in the form 

fi(x) = /„(*) (8) 



where we have defined the functions fi(x) = (x — Xi) 2 + A and f2(x) = ^ X _ x ^ 2+B , and the 
real quantities Xj = Kv gd , A = -Qj = -M n , B = -tt 2 2 = -M 22 , and C = Q 4 C = M 12 M 21 ; 
x here denotes fl The stability profile is determined by the number of real solutions of 
Eq. (jHJ), an integer, say r, between and 4. For absolute stability (for any K, \£jo\), we 
need to have 4 real solutions; in any other case, i.e. if r < 4, the (imaginary part of) the 



6 



4 — r complex solutions determine (s) the growth rate of the instability. Note that x\ 7^ x 2 
expresses the group velocity mismatch v 9y \ 7^ v 9y2 . Negative A (B) means that wave 1 (2) 
alone is stable, and vice versa. 

Let us first consider a wave pair satisfying C > 0, i.e. for M 12 M 2 i ~ P1P2Q12Q21 > 
(a symmetric wave pair, for instance). We shall study the curves representing the functions 
fi(x) and /^(x) on the xy plane. The former one is a parabola, with a minimum at (x\, A). 
The latter one is characterized by a local maximum (for C > 0) at (x 2 , C/B), in addition to 
a horizontal asymptote (the x-axis), since f 2 (x) — > for x — > ±00. Furthermore, for 5 < 
(only), f 2 (x) has two vertical asymptotes (poles) at x = X2 ± v/fsj (see Figs. Now, 
for a stable - stable wave pair (i.e. for A,B < 0), we have seen that the dispersion relation 
(J3J) predicted stability. This result regarded the equal (or vanishing) group velocity case, 
Vg,i — v g,2, and may be visualized by plotting fi(x) and ^(x) for x\ = x 2 and A, B < 0; see 
Fig. Ek- Let us first assume that D = M\\M 22 — M\ 2 M 2 \ = AB — C is positive, implying 
(for B < 0) that A < C/.B; thus, the minimum of fi(x) lies below the local maximum of 
f2(x). Thus, 4 points of intersection exist (cf. Fig. EK), for x\ = x 2 and A,B < 0; this fact 
ensures stability, as we saw above via analytical arguments, for v 9t i = v 9j2 and Mn, M22 > 
(both waves individually stable). Now, considering v g ^ 7^ v g>2 results in a horizontal shift 
between the two curves (cf. Fig. |2t>), which may exactly result in reducing the number of 
intersection points from 4 to 2 (enabling instability). Therefore, a pair of stable waves may 
be destabilized due to a finite difference in group velocity. 

Still for a stable-stable wave pair (A,B < 0), let us assume that D = AB — C < 0, 
implying (for B < 0) that A > C/B. Thus, the minimum of f\(x) here lies above the local 
maximum of f 2 {x), and only 2 points of intersection now exist, (shift the parabola upwards 
in Fig. El to see this); this fact imposes instability (for D < 0), as predicted above. 

Considering an unstable - unstable wave pair (i.e. A > and B > 0) with D = AB—C > 
(A > C/B). Plotting fi(x) and f 2 (x) for X\ = x 2 and A, B > (see Fig. EJ), we see that the 
minimum of fi(x) lies above the local maximum of f2{x). No points of intersection exist, a 
fact which prescribes instability. Considering v 9t \ 7^ v 9j2 simply results in a horizontal shift 
between the two curves, which does not affect this result at all. On the other hand (still 
for A > and B > 0), now assuming that D = AB — C < 0, i.e. A < C/B, results in a 
vertical shift downwards of the parabola in Fig. Et; at least 2 complex solutions obviously 
exist, hence instability. Therefore, a pair of unstable waves is always unstable (V v g ^,v g ^ 2 ). 
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Still for C > 0, one may consider a stable - unstable wave pair (say, for A < and B > 0, 
with no loss of generality): the plot of fi and fa (here omitted) would look like Fig. El upon 
a strong vertical translation of the parabola downwards (so that the minimum lies in the 
lower half-plane, since A < 0). Instability (r = 2) dominates this case also. 

Let us now consider a wave pair satisfying C < 0, i.e. M 12 M 2 i ~ P1P2Q12Q21 = P 2 Q 2 < 
0; this has to be an asymmetric wave pair. Again, different cases may be distinguished. 

For an stable-unstable wave pair, i.e. say for A < < B, different possibilities exist: cf. 
Fig. HJ where 4 points of intersection ensure stability, for v g ^ ~ v g ^)- However, either a 
(horizontal) shift v g difference or (a vertical shift) in A may render the system unstable. 

For A, B > (both waves intrinsically unstable), one easily sees that no intersection 
occurs (figure omitted; simply translate the parabola upwards in Fig. HJ); the pair is unstable. 

Finally, for a stable-stable wave pair (A,B < 0), the wave pair may stable (see Fig. |3J); 
this configuration is nevertheless destabilized either by a velocity misfit (a horizontal shift) 
or a vertical shift (in A). 

In conclusion, we have investigated the occurrence of modulational instability in a pair 
of coupled waves, co-propagating and interacting with one another. Relying on a coupled 
NLS equation model, we have derived a complete set of explicit (in)stability criteria, in 
addition to exact expressions for the critical wavenumber thresholds. Furthermore, we have 
traced the role of the group velocity mismatch on the coupled wave stability. The results 
are readily applied to a set of coupled Gross-Pitaevskii equations (modelling a pair of BECs 
in condensed boson gases), as exposed here, as well as in a variety of physical situations. 
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FIG. 1: The square of the instability growth rate 7 2 = — is depicted versus the perturbation 
wavenumber K (arbitrary parameter values). Notice the difference from the single wave case (lower 
curve) . 
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FIG. 2: The functions f\{x) (parabola) and f2(x) (rational function, two vertical asymptotes) 
defined in the text are depicted, vs. x, for A = B = —1, C = 0.5 (so that D = AB — C = +0.5 > 0), 
xi = X2 = (equal group velocities). Note that a group velocity mismatch (a horizontal shift) may 
destabilize a pair of (stable, separately) waves (i.e. reduce the intersection points from 4 to 2). 
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FIG. 3: The functions fi(x) and /^(a?) are depicted, for A = B = +1, C = 0.5 (so that D = 
AB — C > 0), and xi = X2 = 0. At most 2 intersection points may occur by translation. A pair of 
(unstable, separately) waves is always unstable. 




FIG. 4: Stable-unstable wave interaction: the functions /i(x) and /2(x) are depicted, for A = —1.48, 
B = +1, C = —1.5, and x\ = X2 = 0. This (stable, 4 intersection points) configuration may be 
destabilized either by a horizontal (v g difference) or a vertical (A value) shift. 
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FIG. 5: Stable-stable wave interaction: the functions fi(x) and /2(x) are depicted, for A = — 1, 
5 = —4, C = —1.5 and xi = X2 = 0. This (stable, 4 intersection points) configuration may be 
destabilized either by a horizontal (v g difference) or a vertical (A value) shift. 
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